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Abstract. Results from the modelling of bars in nearly 300 galaxies are used to test pre- 
dictions from theoretical work on the evolution of bars. Correlations are found between bar 
ellipticity and boxiness, between bar strength and normalised size, between the normalised 
sizes of bars and bulges, and between normalised bar size and bulge-to-total ratio. Bars with 
different ellipticities follow parallel lines in the latter two correlations. These correlations 
suggest that, formed with different sizes and ellipticities, bars slow down and grow longer 
and stronger, in agreement with theoretical work. As a consequence, bar pattern speeds 
should become lower with time, and towards galaxies with more prominent bulges. 
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1. Predictions from Theory 

A number of independent results from nu- 
merical simulations suggest that, in the case 
of galaxies with low gas content, bars should 
slow down with time, as a result of angular 
momentum exchange from the inner disc 
to the outer disc or halo. This slowdown 
allows bars to grow in size, capturing stars 
from the disc. In fact, these simulations 
indicate that bars s hould also grow longer and 



strong er (see e.g. Athanassoula & MisiriotisI 
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Berentzen et al. 2006, and references therein). 
An analyti cal treatment of these processes can 
be found in anassou lal (l2003h . Observations 
suggest that the strong bar in NGC 4608 has 
increased in mass by a factor of a:!.?, through 



the capture of a; 13 per cent of disc stars 
(lGadottill2008ah . 

The purpose of this work is to try to verify, 
using imaging data for a sample of nearly 300 
barred galaxies, whether these theoretical pre- 
dictions are fulfilled. The effects of a high gas 
content and other difficulties are briefly dis- 
cussed in Sect.m 



2. l\/leasuring the Properties of Bars 
from Observations 
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In iGadottil (!2008b), the structural properties 
of nearly 1000 massive local galaxies in the 
Sloan Digital Sky Survey (SDSS) were ob- 
tained through careful bar/bul ge/disc decom- 
positiqn, using thejsuDDA code dde Souza et al.l 
120041 lGadotti2008ah . In this sample, there are 
291 barred galaxies, and the code is able to es- 
timate a number of structural parameters, such 
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Fig. 1. Example of a bar/bulge/disc decomposition using budda with one of the galaxies in the 
sample. The top three panels show, from left to right, the original /-band galaxy image, the model 
and residuals. The horizontal line in the galaxy image marks a length of 5 kpc. The residual 
image is obtained subtracting the model from the original image, and is displayed with a narrow 
intensity range, in order to enhance residual sub-structures. In the residual image, darker pixels 
indicate where the galaxy is brighter than the model, whereas whiter pixels indicate where the 
model is brighter than the galaxy. Surface brightness profiles, obtained from cuts along the bar 
major axis, are shown in the bottom panels. The dashed line corresponds to the original image, 
whereas the black solid line corresponds to the total model. Red, blue and green lines refer, 
respectively, to bulge, disc and bar The dotted line in the lower panel shows the residuals (galaxy 
- model). (Taken from ,GadottL2008bi ) 



as bulge effective radius r^, bulge-to-total ra- 
tio B/T, bar effective radius rgbar and length 
(semi-major axis) Leai-, bar ellipticity e, and 
boxiness c. An example of such decomposi- 
tions is shown in Fig.[T] 

The surface brightness profiles of bars are 
modelled with Sersic functions, as are those of 
bulges. However, while bulges have Sersic in- 
dices n typically in the range 1 < n < 6, bars 
have Sersic indices ntar typically in the range 
0.5 < «bar :S 1. In addition, bars are fitted 
as a set of concentric, ge nerahsed ellipses (see 
lAthanassoula et al.l[l9^ : 



where a and b are the bar semi-major and semi- 
minor axes, respectively, and c is the bar boxi- 
ness. If c = 2, then the bar is a perfect ellipse, 
whereas if c > 2 the bar has a more rectangular 
shape. Since bars generally do have rectangu- 
lar shapes, it is important to fit them with c > 2. 
The results presented here were obtained with 
c as a free parameter The bar ellipticity is thus 
1 - b/a. As noted in IGadottil (i2008a). bar el- 
lipticities obtained through detailed image de- 
composition are more reliable than those ob- 
tained via ellipse fits to isophotes. In the latter, 
the superposition of light from bulge and disc 
makes the isophotes in the bar region rounder. 
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Fig. 2. Bar ellipticity plotted against boxiness 
for galaxies with classical and pseudo-bulges. 



leading to a systematic underestimation of the 
bar ellipticity. 

It should also be noted that, due to the rel- 
atively poor spatial resolution of SDSS im- 
ages, bars with Lsai ^2-3 kpc, typically 
seen in very late-type spirals ( later than Sc - 
lElmegreen & ElmegreenlfTQSSl) . are frequently 
missed. Thus, the results presented here con- 
cern bonafide, large bars, as those typically 
seen in early-type disc galaxies. 

3. Bar Strength and Growth 

Figure |2] shows that bar ellipticity and box- 
iness are correlated, and this does not de- 
pend on whether the galaxy hosts a clas- 
sical b ulge or a pseudo-bulge (see iGadottil 
l2008bl for how pseudo-bul ges are identified). 
Furthermore, both parameters are related to 
the strength of bar. Keeping everything else 
the same, a more eccentric (or more rectan- 
gular) bar introduces a more substantial non- 
axisymmetric perturbation in the galaxy poten- 
tial. Thus, one can use the product e x c as a 
measure of bar strength. 

Now, one can see if bar size correlates with 
bar strength (parameterised as e x c), as in the 
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Fig. 3. Normalised measures of bar size plot- 
ted against bar strength for galaxies with clas- 
sical and pseudo-bulges. Longer bars tend to 
be stronger. 



theoretical studies mentioned in Sect. [T] One 
can use either r^ b^ or Lbsi- as a measure of 
bar size. However, larger discs will form larger 
bars, and, in fact, bar size is correlated with 
disc siz e, at least in early-type disc galaxies 
(see e.g. lErwinll2005l) . Therefore, bar size must 
be normalised, and, with this aim, one can di- 
vide r^ bar or Lfiar by the disc scalelength h, or 
by the radius that contains 90 per cent of the 
galaxy light R90, or by the semi-major axis of 
the 24 r-band mag arcsec"^ isophote r24. One 
can obtain h and r24 from the budda models, 
while R90 is available in the SDSS database. 
The results presented here are essentially sim- 
ilar for all such possible parameterisations of 
normalised bar size. Figure[3]shows that longer 
bars tend to be stronger in galaxies with clas- 
sical and pseudo-bulges, in agreement with the 
theoretical work. 

Figure H] shows that the normalised sizes 
of bars and bulges are also correlated. This 
figure confirms the co rrelation found by 
lAthanassoula & Martinet! (Tl 980) with a sample 
of 32 galaxies, and shows that this correlation 
is present in a larger range of normalised bar 
and bulge sizes than previously found. With 
the larger sample of the present study, a new 
aspect of the correlation is found, namely that 
bars with different ellipticities seem to follow 
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Fig. 4. Correlation between the normalised 
sizes of bars and bulges, for bars in three bins 
of ellipticity. The solid line shows a fit to all 
bars and is the same in the three panels. Note 
that when e < 0.6 most of the data points 
are below the line; when 0.6 < e < 0.7 the 
points follow the line more closely; and when 
e > 0.7 most of the points are above the line. 
Although there is no clear correlation for bars 
with e > 0.7 alone, bars with different ellip- 
ticities seem to describe parallel lines in this 
correlation. 



parallel tracks, although there is no clear cor- 
relation for bars with e > 0.7 alone. 

The correlation between the normalised 
sizes of bars and bulges suggests that the 
growth of both components is somehow con- 
nected. Consistent with this idea, Fig.|5]shows 
that the normalised size of bars is also cor- 
related with bulge-to-total ratio. Longer bars 
tend to reside in galaxies with more conspic- 
uous bulges. And, again, bars with different 
ellipticities describe parallel Unes in this re- 
lation. This is now more evident than in the 
relation between the normalised sizes of bars 
and bulges (Fig.H)). Furthermore, it is clear that 
bars of all ellipticities, including those with 
e > 0.7, follow a correlation between their nor- 
malised sizes and BjT . 

Galaxies with larger va lues of BJT ten d to 
be more massive (see e.g. lGadottill2008bh . In 
the current picture of galaxy formation, more 
massive galaxies are believed to form earlier 
than less m a ssive galaxies, as suggested by 
ICowie et al.l (Il996i) . In this case, BjT may 
serve as a proxy for time: galaxies with larger 
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Fig. 5. Normalised bar size plotted against 
bulge-to-total ratio, for bars in three bins of el- 
lipticity. The top panel shows the data points 
and the solid line is a fit to all points. Longer 
bars tend to reside in galaxies with more con- 
spicuous bulges. Furthermore, bars with differ- 
ent ellipticities describe parallel lines in this re- 
lation. This is better seen in the lower panel, 
which shows separate fits in the different ellip- 
ticity bins, with the slope fixed at the value in 
the fit to all points. 



BjT are older than those with smaller BjT . 
If this is correct, then the correlation between 
normalised bar size and BjT (Fig.|5j indicates 
that bars are longer in more evolved galaxies, 
and thus grow longer with time. Combining 
this result with the correlation between nor- 
malised bar size and bar strength (Fig.O, one 
concludes that bars grow longer and stronger 
with time, in agreement with the theoretical 
predictions described in Sect.[T] 

Us ing a diffe rent methodology, 
Elmegr een et al.l (|2007|) also find that nor- 
malised bar size correlates with bar strength 
and with the galaxy central densi t y, and reach 
similar conclusions. ISheth et alJ (l2008h find 
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that more massive galaxies have their bars in 
place at higher redshifts, whereas less massive 
galaxies form bars at later times. They suggest 
that discs in more massive systems reach a 
dynamical maturity earlier than those in less 
massive systems, and thus are able to form 
bars at earlier times, in agreement with the 
results and interpretations given here. 

The current understanding of the orbital 
structure in barred galaxies, and the obser- 
vation that most bars en d near their coro- 
tation radius (see e.g. Rau tiainen et al.l l2005l 
and references therein), tell us that bars gen- 
erally cannot grow longer if they do not 
slow down. Therefore, the results presented 
here also suggest that bars slow down with 
time, which is, again, consistent with the- 
ory. Such a relation between the bar pattern 
speed Q and BjT is in the sense that bars 
rotate slower in galaxies with more promi- 
nent bulges, since these galaxies have longer 
bars. This relation can also be seen as a 
dependence of Q. on Hubble type, although 
there is some scatter in the relation between 
^/r and Hubble type (e.g. L aurikainen et al.l 
l2007HGraham & Worlevll2008h . Given that di- 
rect measurements of bar pattern speed are dif- 
ficult, particularly for late-type spirals, it is 
not surprising that no solid conclusions can 
presently be drawn about a dependence of bar 
pattern speed on Hubble type, on direct ob- 
servational grounds ( see lOerssen et al.l l2003t 
lTreuthardtetani2007l) . 

An oustanding and unforeseen new result 
in Figs. |4] and |5] is the existence of parallel 
tracks in e.g. the correlation between bar nor- 
malised size and BjT, for bars with different 
ellipticities. A straightforward way of inter- 
preting the existence of these parallel tracks is 
to conceive that bars should form with different 
normalised sizes and ellipticities, and then fol- 
low a somewhat parallel growth. This is a new 
aspect that can be investigated with theoretical 
work. For instance, can simulations form bars 
with various normalised sizes and ellipticities? 
Do simulations show that, although formed 
with somewhat different properties, such bars 
follow a similar evolutionary path? 



4. Discussion 

The correlations presented above essentially 
corroborate the picture provided by theoretical 
work on the formation and evolution of bars, 
at least when the galaxy gas content has lit- 
tle eff'ect on the bar. The role of gas is still 
a matter of debate (se e Bourn aud et al.|l2005i: 
iDebattista et"al] 120061: .Berentzen et al.ll2007l) . 
Nevertheless, most simulations show that the 
efifect of gas is to weaken the bar with time. 
Furthermore, in some gas-rich simulations, the 
bar strength oscillates substantially, alternating 
between very weak and very strong in cycles. 
The evolution of the bar pattern speed seems 
not to be significantly affected by the presence 
of gas. Bars still slow down in gas-rich simula- 
tions, which means they can also grow longer. 

Another issue that can complicate the inter- 
pretation of the observational results presented 
here is related to the bar vertical buckling in- 
stability. Many simulations show that, at early 
stages, after becoming very strong, bars grow 
vertically thick in their central parts, ofi" the 
disc plane. Due to this process, bars weaken 
dramatically, and then start growing stronger 
again. The strongest buckling instability occurs 
soon after the formation of the bar, and ends 
when the bar is 2 - 3 Gyr old. 

It should also be mentioned that, depend- 
ing on the initial conditions, the evolution of 
bar properties in simulations can be very slow. 
All these issues might contribute to the scat- 
ter seen in the observed correlations above. 
Nonetheless, since these correlations are based 
on the observation of typically large bars, and 
since short bars are associated with gas-rich 
systems, it is likely that most of the bars in the 
sample have not been significantly aff'ected by 
the effects induced by gas. Furthermore, since 
most long bars are not expected to have been 
recently forme d (this comes not only from sim- 
ulations - see Gadott i & de Souzal l2006'). it is 
also likely that most bars in the sample have 
already gone through the most significant ver- 
tical buckling. If these suppositions are correct, 
then one indeed expects to see the correlations 
presented. A similar work, with a sample of 
short bars, in gas-rich systems, is likely to shed 
light on these issues. 
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Finally, it should be noted that the deter- 
mination of the bar parameters involved in the 
results above, namely bar length, ellipticity 
and boxiness, is difficult. The scatter observed 
might at least partially result from measure- 
ment uncertainties (which are themselves also 
difficult to estimate). 
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